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Abstract. We use computer algebra to determine the ranks of arrays of size 
2x2x2 and 2x2x2x2 with entries in the set {0, 1} regarded as a field 
with two elements, as a Boolean algebra, and as non-negative integers. In the 
field case we also determine the canonical forms of the arrays with respect to 
the action of the direct product of the general linear groups. 



1. Introduction 

Multidimensional arrays and the related topic of hyperdeterminants have con- 
nections with algebraic geometry and representation theory, and applications in 
numerical analysis, signal processing, chemometrics and psychometrics. For the 
connections with algebraic geometry see [5j \6[ [10] ; for a connection with represen- 
tation theory see [T]. For surveys of the applications, see pi 151 1^ [TT]. 

Most research on this topic assumes that the arrays have entries in K or C, the 
fields of real and complex numbers. In this paper, we consider arrays with entries 
in {0, 1}, which can be regarded as the field F2 with two elements (1 -|- 1 = 0), as a 
Boolean algebra (1 -f 1 = 1), or as non-negative integers (1 -|- 1 = 2). In this case, 
the total number of arrays is finite, and the problem of determining the rank of an 
array reduces to combinatorial enumeration which can be performed by computer. 
We obtain a classification by rank of all 2 x 2 x 2 and 2x2x2x2 arrays in these 
three cases. In the field case we determine the canonical forms of the arrays with 
respect to the action of the finite groups GL2(F2)'^ and GL2(F2)^ and the extended 
groups GL2(F2)3 x ^3 and G'i2(F2)'^ x ^4. 

The canonical forms of 2 x 2 x 2 arrays have been determined over the real 
numbers [3] and over the complex numbers [H [5] . Analogous results for 2 x 2 x 2 x 2 
arrays over R or C have not yet been found, but see [7] . Since the results for 2x2x2 
arrays over R and C are similar to the results in this paper over F2, we hope that 
our results for 2x2x2x2 arrays over F2 will provide some useful information 
towards a classification of canonical forms of 2x2x2x2 arrays over R and C. 

2. Preliminaries 

We consider n-dimensional arrays of size 2 x • ■ • x 2 (n factors, n = 3,4) with 
entries in {0, 1}. 

Definition 1. The flattening of the array X — 1 < «i, . . . , i„ < 2, is 

fiat(X) = [ a;i...i ••■ Xi^.-.i^ ••• a;2...2 ], 
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where the entries are in lexicographical order of the n-tuples of subscripts: ii ■ ■ ■ in 
precedes i'l ■ ■ ■ i'^ if and only if ij < i'j where j is the least index with ij ^ i'^ . 

Definition 2. If X and Y are two arrays then X precedes Y if flat(X) precedes 
flat(y) in lexicographical order: that is, 2;^^...^^ < yi-^...i^ where ii ■ • • i„ is the least 
n-tuple with Xi-^...i^ ^ yi^...i^. The minimal element of a set of arrays is defined 
with respect to this total order. 

Definition 3. There are three nonzero 2-dimensional vectors: [0,1], [1,0], [1,1]. 
The outer product X = Vi® ■ ■ ■ ®Vn of n vectors Vj — [uji, Wj2], 1 < i < is 

Definition 4. The rank oi X — {xi^...i^) is the minimal number R of terms in the 
expression of X as a sum of outer products: 

The definition of addition depends on the algebraic structure of {0, 1}. For the field 
with two elements, 1 + 1 = 0; for the Boolean algebra, 1 + 1 = 1; for non- negative 
integers, 1 + 1 = 2, and this means that we exclude any sums in which two outer 
products both have an entry 1 in the same position. 

Lemma 5. The only array of rank is the zero array. An array of rank 1 is the 
same as an outer product of nonzero vectors, and there are 3" such arrays. The 
total number of n- dimensional 2 x ■ • • x 2 arrays with entries in {0, 1} is 2^ . 

Algorithm 6. Fix a dimension n. Assume that we have already computed the 
arrays of rank r. To compute the arrays of rank r + 1, we consider all sums X + Y 
where rank(A) = r and rank(y) = 1. Clearly rank(A + Y) < r + 1, but it is 
possible that rank(A + Y) < r, so we only retain those X + Y which have not 
already been computed: the arrays which have rank exactly r + 1. This algorithm 
is presented in pseudocode for rt = 3 in Table [T] 

If {0, 1} is the field F2 with two elements then we consider the action of the 
direct product of n copies of the general linear group GL2(F2). 

Definition 7. The group GL2{¥2) consists of six matrices in lexicographical order: 
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1 1 
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Lemma 8. The group GL2{^2) is isomorphic to S3 permuting the nonzero vectors. 

Definition 9. The small symmetry group of 2 x • • • x 2 arrays over F2 is 
the direct product Gi2(F2)" acting by simultaneous changes of basis along the n 
directions. The large symmetry group of these arrays is the semi-direct product 
0/72(1^2)" XI Sn where Sn acts by permuting the n copies of GL2(F2). The element 
A S GL2{¥2) acts along the first direction of an array X = (xjj...^^) as follows: for 
each (n— l)-tuple 12 ■ ■ ■ in we consider the column vector 

^2---jii ~ [xii2---inj X2i2---in] & ^2^ 

and compute AVi^-.-i^ — [yii2 - j„ 1 2/2i2 - jn]* € ^2; then we define the array A - X io 
heY— The actions along the other n — I directions are similar. 
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f latten(a;) 

return([a;iii, Xii2, Xi2i, X122, X211, X212, X221, 2:222]) 

outerproduct(a, b, c) 

for i = 1, 2 do for j = 1, 2 do for fc = 1, 2 do: 

set Xijk ^ ttibjCk 
return( x ) 

• set vectors ^ { [1, 0], [0, 1], [1, 1] } 

• set arrayset[0] ^ {[0,0,0,0,0,0,0,0]} 

• set arrayset[l] <— { } 

• for a in vectors do for b in vectors do for c in vectors do 

— set X <— f latten(outerproduct(a, 6, c)) 

— if X ^ arrayset[0] and x ^ arrayset[l] then 

set arrayset[l] arrayset[l] U {x} 

• set r ^ 1 

• while arrayset[r] 7^ { } do: 

— set arrayset[r+l] { } 

— for X £ arrayset[r] do for y £ arrayset[l] do 

* set z ^ [xi+yi, . . . , xs+ys] 

* if z ^ arrayset[s] for s = 0, . . . , r+1 then 

set arrayset[r+l] arrayset[r+l] U {z} 

— set r ^ r + 1 

• set maximumrank ^ r — 1 

Table 1 . Algorithm [6] in pseudocode 



Lemma 10. The actions of the symmetry groups do not change the rank. 

Proof, de Silva and Lim [Sj Lemma 2.3, page 1092] applies to any field. □ 

The actions of the symmetry groups decompose the set of n-dimensional arrays 
into a disjoint union of orbits; the arrays in each orbit are equivalent under the 
group action. 

Algorithm 11. Fix a dimension n and consider 2 x • • • x 2 arrays over F2. Assume 
we have computed the set of arrays for each possible rank, and that these sets are 
totally ordered. For each rank, we perform the following iteration: 

• Choose the minimal element of the set of arrays. 

• Compute the orbit of this element under the action of the symmetry group. 

• Remove the elements of this orbit from the set of arrays of the given rank. 

This iteration terminates when there are no more arrays of the given rank. This 
algorithm is presented in pseudocode for n = 3 in Table [5] 

Definition 12. The minimal element in each orbit is the canonical form of the 

arrays in that orbit. 
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unf latten(a;) 

set t 

for i = 1, 2 do for j = 1, 2 do for fc = 1, 2 do: 

set t <— t + 1 

set Uijk <- xt 
returii(2/) 

groupaction((7, x, m) 

set y <r- unf latteii(x) 

if m = 1 then 

for j = 1,2 do for fc = 1,2 do 
set V ^ [yijk, y2jk] 
set w [giivi + 512W2, g2ivi + 322^2 ] 
for i = 1, 2 do: set yijk ^ Wi 

if m = 2 then . . . (similar for second subscript) 

if TO = 3 then . . . (similar for third subscript) 

return( f latten( y ) ) 

sinallorbit(a;) 

set result ^ { } 
for a e GL2(F2) do: 

set y groupaction(a, x, 1) 
for b e GL2{¥2) do: 

set z •(— groupaction(&, y, 2) 
for c G GL2(F2) do: 

set w <r- groupaction(c, z, 3) 
set result ^ result U {w} 
return( result ) 

largeorbit(x) 

set y <r~ unf latten(a;) 
set result { } 
for p £ S3 do: 

for i = 1, 2 do for j = 1, 2 do for fc = 1, 2 do: 
set m <— [i, j, k] 

set Zijk 2/mp(i)mp(2)mp(3) 

set result <~ result U sinallorbit(f latten( 
return( result ) 

for r = 0, . . . , maximumrank do: 

set representatives[r] <— { } 

set remaining <— arrayset[r] 

while remaining ^ { } do: 
set X <— remaining[l] 
set xorbit <r- largeorbit(a:) 
append xorbit[l] to representatives[r] 
set remaining <— remaining \ xorbit 

Table 2 . Algorithm [11] in pseudocode 
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rank orbit size canonical form 
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■ 



■ 






1 27 


■ 



■ 

1 




2 54 
2 108 
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■ 
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1 ■ 

1 

1 ■ 






3 54 
3 12 


■ 

1 

■ 1 

1 


1 ■ 

1 

1 ■ 
1 1 





Table 3. Large orbits of 2 x 2 x 2 arrays over F2 



Lemma 13. Lower hounds for the number of canonical forms for the small sym- 
metry group and the large symmetry group are respectively 



-22"- 




- 22" - 


6" 




6"n! 



Proof. Lemma [5] shows that there are 2^" such arrays. Definition IHlimphes that the 
small and large symmetry groups have orders 6" and 6"n! respectively. The claim 
follows from the theory of group actions on a finite set. □ 

Remark 14. For 3 < n < 6, we have the following lower bounds for the number of 
orbits for the small and large symmetry groups. From this it is clear that complete 
results will not be publishable for n > 5: 

n lower bound (small group) lower bound (large group) 

3 2 1 

4 51 3 

5 552337 4603 

6 395377745064077 549135757034 



3. Arrays of size 2x2x2 

The set of 2 x 2 x 2 arrays with entries in {0, 1} contains 256 elements. We 
represent such an array X — {xijk ) in the matrix form 

a;ii2 2^122 

a;212 2^222 ' 



Mat(X) = 



a;iii X121 
X211 X221 



where the third subscript distinguishes the left and right blocks, which are the first 
and second frontal slices. 
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rank ones number representative 
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4 4 
4 5 



Table 4. Ranks and minimal representatives for 2 x 2 x 2 Boolean arrays 



3.1. The field with two elements (1 + 1 = 0). Algorithm [6] shows that in this 
case the maximum rank is 3; the number of arrays of each rank is 1, 27, 162, 66. 
The percentages of ranks 0, 1, 2, 3 are approximately 0, 11, 63, 26; in contrast, the 
percentages over M are approximately 0, 0, 79, 21. For the large symmetry group 
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GL2(F2)'^ X 5*3, the ranks, orbit sizes, and canonical forms are given in Table [S 
For the small symmetry group GL2{¥2)^, the first orbit in rank 2 splits into three 
orbits each of size 18 with canonical forms 
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For the small symmetry group, there are eight orbits, the same as in the real case. 

3.2. The Boolean algebra (1 + 1 = 1). The maximum rank is 4; the number 
of arrays of each rank is 1, 27, 130, 88, 10. The percentages of ranks 0, 1, 2, 3, 4 
are approximately 0, 11, 51, 34, 4. Instead of canonical forms for a group action, 
which do not exist in the Boolean case, we partition the arrays in each rank by the 
number of entries equal to 1; the results are given in Tabled 

3.3. Non-negative integers (1 + 1 = 2). The results are the same as in the 
Boolean case. This has the corollary that every 2x2x2 Boolean array of rank 
r can be written as the sum of r outer products such that no two terms have an 
entry 1 in the same position. (When we consider arrays of size 2x2x2x2, the 
Boolean and integer cases will no longer be identical.) 



4. Arrays of size 2x2x2x2 
The set of 2x2x2x2 arrays with entries in {0, 1} contains 65536 elements. 

4.1. The field with two elements (1 + 1 = 0). Algorithm [5] shows that in 
this case the maximum rank is 6. The number of arrays of each rank and the 
approximate percentages are as follows: 

rank 1 2 3 4 5 6 

number 1 81 2268 21744 37530 3888 24 

0.002 0.124 3.461 33.179 57.266 5.933 0.037 

For the large symmetry group GL2{¥2)'^ x S4, there are 30 orbits; the ranks, orbit 
sizes, and canonical forms are given in Table [5l For the small symmetry group 
GL2(]F2)'*, there are 112 orbits. The large orbits split into small orbits as follows, 
where we mention only those large orbits that are not small orbits, and write 
X ^ y ■ z to indicate that large orbit x splits into y small orbits each of size z: 
rank 2 rank 3 rank 4 rank 5 



3^6 


54 


6^4- 


162 


15 - 


■>4- 


648 


26 6 • 108 


4^4 


324 


7^4- 


36 


16 - 


->4- 


1296 








8^6- 


648 


17- 


^ 3- 


36 








9-^4- 


648 


18 - 


-> 3- 


324 








10 ^ 4 


•648 


19 - 


^6- 


324 








11 ^ 3 


• 1296 


20 - 


^ 3- 


648 








12 ^ 6 


• 1296 


21 - 


^ 12 


•648 












22 - 


^ 6- 


216 












23 - 


^6- 


1296 












24 - 


^3- 


1296 












25 - 


^6- 


648 





8 



BREMNER AND STAVROU 



4.2. The Boolean algebra (1 + 1 = 1). The maximum rank is 8. The number of 
arrays of each rank and the approximate percentages are as fohows: 

rank 012 3 4 5678 

number 1 81 1804 13472 28904 17032 3704 512 26 

«% 0.002 0.124 2.753 20.557 44.104 25.989 5.652 0.781 0.04 

The results are given in Table [6] by rank and number of entries equal to 1. 

4.3. Non-negative integers (1 + 1 = 2). The maximum rank is 8. The number 
of arrays of each rank and the approximate percentages are as follows: 

rank 012 3 4 5678 

number 1 81 1756 12848 28788 17568 3908 560 26 

w % 0.002 0.124 2.679 19.604 43.927 26.807 5.963 0.854 0.04 

The results are given in Table [7] by rank and number of entries equal to 1. 
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24 
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3888 
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5 


648 
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27 
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648 


0001011001101000 


28 
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1296 


0001011001101001 


29 
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24 


110 1 1 1 10 11110 1 



Table 5. Large orbits of 2x2x2x2 arrays over F2 
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26 
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8 
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9 
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35 
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36 


4 
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12 
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13 
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4 


14 
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16 
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6 
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54 


6 


10 


848 
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160 
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9 


16 


0110100110010111 
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8 


10 


8 


0110101111010110 



Table 6. Minimal representatives for 2 x 2 x 2 x 2 Boolean arrays 
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21 


3 


8 


1664 


0000000111101111 


22 


3 


9 


864 


0000001111011111 


23 


3 


10 


608 


0000011011111111 


24 


3 


11 


384 


0000011111111111 


25 


3 


12 


256 


000 1101111111111 


26 


3 


13 


96 


000 1111111111111 


27 


3 


14 


32 


11111111111111 


28 


4 


4 


228 


0000000001101001 


29 


4 


5 


1648 


OUUUUUUUUilUlUli 


30 


4 


6 


4048 


00000001001 11 110 


31 


4 


7 


5856 


U U () U U U i U i i U i 1 i i 


32 


4 


8 


6544 


(J 1 1 1 1 1 1 1 1 


33 


4 


9 


5104 


0000001101111111 


34 


4 


10 


3056 


0000011110111111 


36 


4 


11 


1504 


0000111111110111 


36 


4 


12 


448 


0001011111111111 


37 


4 


13 


256 


0011110111111111 


38 


4 


14 


80 


110111111111111 


39 


4 


15 


16 


111111111111111 


40 


6 


5 


128 


0000000110010110 


41 


6 


6 


1008 


000000011001011 1 


42 


6 


7 


2416 


0000001101101101 


43 


5 


8 


3568 


0000011001101111 


44 


5 


9 


4304 


(0 1 1 (H) 1 1 1 1 1 1 1 


45 


5 


10 


3088 


1110 1111111 


46 


5 


11 


1984 


10 1111110 1111 


47 


5 


12 


904 


(J 1111111110111 


48 




13 


160 


Ui--LU-LUi--Li-i--l--L-Li--l--L 


49 


6 


14 


8 


0111111111111110 


50 


6 


6 


56 


0000011001101001 


51 


Q 


7 


448 


110 110 10 11 


52 


6 


8 


848 


1110 11110 1 


53 


6 


9 


928 


10 110 110 1111 


54 


6 


10 


1040 


000 10 110 01111111 


66 


6 


11 


368 


0001011101111111 


66 


6 


12 


172 


0110101111011111 


57 


6 


13 


48 


0110111111110111 


58 


7 


7 


16 


10 1 1 1 10 1 1 


59 


7 


8 


128 


10 1 1 1 10 1 1 1 


60 


7 


9 


160 


10 1111110 10 1 


61 


7 


10 


112 


0011110111010110 


62 


7 


11 


128 


0110100110111111 


63 


7 


12 


16 


0110101110111111 


64 


8 


8 


2 


0110100110010110 


66 


8 


9 


16 


0110100110010111 


66 


8 


10 


8 


0110101111010110 



Table 7. Minimal representatives for 2x2x2x2 integer arrays 



